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An investigation has been conducted to' determine' the effects of 
j et-to-f ree-sf ream dynamic pressure ratio on the penetration and mixing 
of a sonic hydrogen jet injected normal to a Mach 4 airstream. The 
hydrogen gas was injected from a circular nozzle flush mounted in a 
flat plate with, a turbulent boundary-layer thickness of 2. JO injector 
exit diameters at the injector station. The investigation was conducted 
for values of the dynamic pressure ratio ranging from 0.5 to 1.5. At 
five downstream stations between J and 200 injector diameters the mixing 
region was surveyed to obtain hydrogen volume concentration and pressure 
profiles. 

Results of the investigation indicate that' the thick boundary 
layer had significant effects on the penetration and maximum concentra- 
tion trajectories when compared to data correlations from other sources. 
The penetration trajectory was found to be propcprtional to the dynamic 
pressure ratio raised to the 0.3 power. The decay of the maximum con- 
centration was very rapid, in the near field and inversely proportional 
to (x/dj ) * at downstream distances greater than 30 jet diameters. 
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At a particular value of x/d the maximum concentration 'was propor- 
tional to q, Uondimensiona! concentration profiles, represented 

hy Gaussian-type functions, on the vertical centerline showed similarity 
at values of x/d equal to or greater than 60. 
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v. introduction. 

Recent projects of the National Aeronautics and Space Administra- 
tion have "been concerned with the design and development of research 
engines for operation in the Mach 4 to -Mach ■ 8 flight regime (ref. l). 

The engine is an axisymmetric hydrogen-fueled ramjet and employs' super- 
sonic combustion at flight speeds above Mach 6. Advanced hypersonic 
vehicles, such as a reusable launch. vehicle with an airbreathing pro- 
pulsion system for the first stage., are currently under investigation 
(ref. 2)-. The propulsion system will he a supersonic combustion ramjet 
with operation up to Mach 12. A, flight Mach number of 12 corresponds 
to a combustor-entrance Mach number of h. Preli m inary designs indicate 
that with the. length of the inlet and compression surfaces required at 
these high flight Mach numbers, a significant portion of the flow, 
entering the combustor will consist of a boundary layer. The accurate 
analysis and design of a supersonic combustor requires a knowledge of 
the fuel-air mixing characteristics. In addition, injecting the fuel in 
a way that provides for a nearly uniform fuel distribution and a short 
mixing length without producing significant thrust penalties, is 
desirable. Parallel or coaxial injection, while contributing signifi- 
. cantly to the thrust, requires a longer mixing length and has less 
penetration than does injection perpendicular to the airstream (refs.. 3 
and V). 

Normal sonic injection of various gases and gaseous mixtures 
issuing from a discrete circular hole in a flat plate into a supersonic 
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airstream has been investigated and reported in references 4 to 9 
inclusive. These data are generally for conditions corresponding to 
high values of j et-to-f r,ee-stream dynamic pressure ratio. Unpublished 
data available at the Langley Research Center indicate that the short- 
est mixing length is obtained at the lowest value of j'et-to-free- stream 
dynamic pressure ratio at -which the injector operates choked. Analyti- 
cal methods for predicting the initial penetration of the jet and the 
jet trajectory in the unconfined supersonic main-stream have been 
developed from empirical or semi— empirical data correlations (see 
refs. 7 to 9)* Generally, ‘ any effect of the molecular 'weight of the 
injected gas or the main-stream boundary-layer thickness has not been 
considered. For application to the design of a supersonic combustor, 
the primary interest is in the .far-field mixing region rather than the 
complex flow in the vicinity of the injector. 

The present investigation was conducted to provide information 
about the effect of jet-to-free-stream dynamic pressure ratio and' a 
thick- boundary layer on the mixing of hydrogen injected normal to a 

* t 

uniform supersonic airstream. * These data are needed to aid in con- 
structing analytical methods that yield reliable predictions about the 
mixing process in supersonic combustors . Hydrogen gas was injected 
from a 0.10l6-centimeter-d±ameter sonic nozzle perpendicular to the 
surface ( of a flat plate mounted in' the 23-centimeter-square tunnel test 
section. " The tests- were conducted at a free-stream Mach number of 4.03, 
stagnation temperature of 300° stagnation pressures of 13.6 and 
20.4 atmospheres, giving Reynolds numbers per meter of 6.19 X 10? and 
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9*28 X 107, respectively. Boundary-layer thickness on the flat plate 
at the injector station was approximately 2.70 injector diameters. Jet 
conditions were such that ratios of jet-to-free-stream momentum flux of 
0.5, O.75, 1.0, and 1.5 were obtained. Measurements of hydrogen volume 
concentratipn} pitot pressure, and static pressure were obtained by 
vertical and horizontal surveys of the flow field at downstream stations 

i 

of 7, 30, 60, 120, and 200 injector diameters. 



VI. LIST OF SYMBOLS 


streamiube cress-section area, meter^ 
exponent defined in equation (9) 
exponent defined in equation (10) 
drag coefficient. 

injector nozzle exit diameter, meter 
equivalent Jet exit diameter = dK 1 / 2 
fuel -air mass ratio 
injector nozzle discharge' coefficient 
mass flow rate, kilogram/second 
Mach number 

absolute' pressure, Newton /met er2 or atm 
effective back pressure, Newton/meter^ 

ratio of Jet-to-free-stream dynamic pressure = .(pv 2 )^/(pv 2 ) 00 
Reynolds number 

distance measured along centerline of emerging Jet, meter 
absolute, temperature, °Kelvin 
velocity, meter-/ second 

longitudinal coordinate 

» x 

lateral coordinate 
vertical coordinate 

» > i 

hydrogen mass fraction 

air mass flux parameter = (pv) x (l - a.) kilogram/met er^- second 
boundary -layer thickness, meter 
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'5* '■ boundaiy -layer .displacement thickness, meter 

v t " hydrogen volume fraction 

1 * i \ 

' i , hydrogen mass flow rate parameter = a(pV ) x / (pV )j 

* , * 

9 ’slope of emerging get centerline measured from 

horizontal^ degree 

.p : mass density, kilogram/meter3 

f 

Subscripts: 

* ? 

00 free-stream conditions 

0 , conditions at the edge of the mixing region where v = 0.005 

1 conditions in undisturbed flow upstream -of injector 

5 • conditions at which the mass concentration is half maximum 

j jet conditions 

t stagnation conditions 

x survey point 

max maximum value 

ref .reference value 

a conditions at which concentration is maximum - 

■Superscripts: 

) average . quantity 

( )’ conditions behind normal shock 



VII. ' MODEL AND FACILITY 


Test Apparatus and Model 

A sketch of the stainless' steel rectangular flat plate used in 
the experiments is shown in figure 1. The plate leading edge was a 
2° wedge tapering to a cylindrical leading edge of approximately 
0.0127 centimeter thickness, followed "by a 10° wedge. A 0.1016- 
centimeter-diameter sonic nozzle was flush mounted perpendicular to 
the- plate surface 18.6 centimeters from the. plate leading edge. Details 
of the nozzle are shown inset in figure 1. The injector tube had a 
constant area section approximately three exit diameters long and was 
fed by a section of tube .4- centimeters long and 0.1524 centimeters 
inside diameter that acted as a plenum. * The plenum section was fed by 
a 0.476 -centimeter -diameter tube in which the jet total pressure was 
measured. For sonic operation of the injector, the Mach number in the 
plenum section is approximately O.3O .and would probably be fully 
developed pipe flow. Because of the rapid acceleration, the flow' in 
the injector tube would be only .slightly affected by the boundary layer 
in the supply tube and would not be fully developed pipe flow.. 

The tests were conducted in a continuous flow supersonic tunnel 
with the flat plate spanning the 25-centimeter by 25-centimeter test 
section. The top surface of the plate was positioned horizontally ' 
along' the longitudinal centerline of the tunnel test section. The 
tunnel - exhausted to atmosphere 1 and had a two- dim ensional fixed geometry 
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nozzle and a second minimum followed by a subsonic diffuser. The test 
section Mach number above the boundary layer at the injection station 
was it-. 03 (ref. it-).' Tests were conducted at stagnation temperatures of 
300° K, stagnation pressures of 13*6 and 20.4 -atmospheres , giving 
corresponding unit Reynolds numbers of 6. 19 X 107 and 9.28 X 10? per 
meter. Boundary-layer surveys were made at the injector exit station 
at both of the tunnel stagnation pressures. Profiles of total pressure 
and velocity for both cases are presented- in figure 2. The boundary- 
layer thickness was taken 'as the point; at which the velocity reached 
99 percent of free stream with values of 2.76 injector exit diameters 
obtained at both free-stream stagnation pressures. Also presented in 
figure 2 are theoretical results obtained from a computer program 
reported in reference 10. The theoretical results were obtained for 
conditions corresponding to a free-stream total pressure of 20.4 atmos- 
pheres and agreed well with- the data. 

Secondary flow . - The routing of the hydrogen gas within the test 
cell is presented schematically in figure 3. Hydrogen gas was stored 
in trailers and supplied to the test cell by a 7* 62-centimeter-diameter 
line equipped with pressure-controlled valves and a nitrogen purge 
system. An electrically controlled, air-operated three-way ball valve 
was used to shut off hydrogen flow inside the test cell and to vent the 
supply line within the building to the atmosphere. Hydrogen flow rate 
was measured by a 0.3175“ cen 4imeter-diameter sharp-edged comer tap 
orifice meter. Upstream pressure at the orifice meter was set and 
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maintained by an air-controXXed pressure-regulated valve and the 
orifice pressure drop and hydrogen jet total pressure controlled by an 
air-operated valve immediately downstream of the flow meter. Hydrogen 
total temperature was measured with a standard iron- const ant an thermo- 
couple inserted in a filter positioned between the -orifice meter and - 
injector. Jet total pressure was. measured near the injector exit by- 
mounting a 0. 1016-centimeter out side -diameter tube in the- O.476 
centimeter-diameter injector supply tube. .Calculations indicated that 
the measured pressure would be within 99 percent of' the jet total 
pressure. The apparatus was operated over a jet total pressure range 
of 2 to 4 'atmospheres corresponding to jet-to -free- stream dynamic 
pressure ratios of O'. 5 to 1.5. The exact test conditions are presented 
in the following table. 


Test 

condition 

Q.r 


p t,j , 
(atm) 

% 2 

(gm/sec ) 

‘ Re '* 

- (m" 1 ) ‘ 

5 /d 

1 

0.50 

' 20 . 4 ' 

1.986 

0.0820 . 

9.28 x 10T 

' 2.70 

2 

0.75 

20.4 

‘2.980 

0 

H 

'OJ 

0 

r t 'i 

9.28'x io 7 

2.7O 

*3 

1.00 

13.6 

i- 

■ 2 . 647 - 

O.IO94 

6.19 x IQ 7 ' 

2. 70 

4 

1.50 

13.6 

3 '. 96o 

0.1641 

6.19 x io 7 * 

2.70 


Chromatograph system . - A tube attached" to the inj ector nozzle 
supply line supplied 100 percent hydrogen .samples for full-scale 
chromatograph readings. The sample collection , and analysis system is 
shown in figure 3* During a survey gas samples of the hydrogen-air 
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mixture were taken, through a pitot probe with the aid of a vacuum pump 

J 1 

at mass flow rates as high as 4-17 x 10“^ gm/sec. She sample flow to 
the chromatograph was metered to 5*56 x 10"^ gm/sec hy an electrically- 
controlled microvalve and the remaining flow bypassed and exhausted 
into the- test cell. Both the sample and bypassed mass flow rates were 
measured by thermo conductivity flow meters. A nitrogen bottle was used 
to insure complete purging of the gas collection and analysis system. 

Instrumentation 

G-as analyzer .- The volumetric concentration of hydrogen in the gas 
samples was measured by a process gas chromatograph. The sample gas 
■and a carrier gas (nitrogen) flow continuously through the chromato- 
graph. - At the beginning of a 1-minute cycle a portion of the sample 
gas is isolated and forced by the carrier gas through a molecular sieve 
and a column consisting of a length of stainless steel tubing packed 
with silica .gel. This provides a qualitative identification of each 
component since each will process through the column at a predictable 
rate. The quantity of each component is determined by four thermo- 
conductivity detectors of which two are always exposed to the carrier 
gas. The unbalance of the detector bridge provides a voltage output 
proportional to the cooling effect, and hence, is a measure of the 
concentration of the separated sample- components relative to the 
carrier gas. The voltage, output is Recorded by a pen deflection on 
a strip chart. Readout controls were -adjusted so' that only the hydrogen 
concentration was detected. The 'pen- deflection for 100 percent hydrogen 
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taken from the supply line -was recorded before each’ survey, and the 
repeatability of the instrument checked to a variation of less than 
l/2 percent full scale. This corresponds to an error in hydrogen 
volume fraction of ±0. 005. Further information about gas chromatography 
may be found' in references 11' and 12.. The gas analyzer was calibrated 
■with, known mix tures of hydrogen and air. The calibration points and 
an equation for a curve through the points are presented in figure 4. 
During the course of the tests, the calibration was kept up to date by 
spot-check calibrations, using hydrogen -nitrogen mixtures of known 
concentration. 

Probe description . - The gas-sampling pitot probe and the' static 

probe are shown in figure 5* The pitot -sampling probe is a boundary- 

layer survey type with the probe tip mounted in' a 7 • 9^- rrnl Tlm ef.er- 

diameter supporting tube offset to -allow for actuator rod clearance.. 

The actuator mechanism provided for probe movement for vertical', 

traversing and yaw in the horizontal; plane. A variable resistance pot 

electrically coupled to a counter indicated the probe position. The 

' * 

counter was calibrated with a precision dial gauge and gave the- probe 
position with an accuracy of ±0.127 millimeter in the traverse mode 
and ±0.10 in the yaw mode. The static pressure probe was of similar 
design and had a 28° cone angle and four 0. 203 -millimeter orifices 
located at 14 probe' diameters from the tip. 

Flow measurement . - The injected gas was measured with a 0.3175 
centimeter-diameter sharp-edged comer tap orifice meter as shown in 
figure 3* Orifice meter upstream static pressure and pressure drop 
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were measured "by a 1.379 X 10^ Ij/m 2 transducer and a 6.895 x lO^ N/m 2 

' 1 , t 

differential pressure transducer, respectively. The static temperature 
at the meter was assumed to be the same as the jet stagnation tempera- 
ture. Hydrogen mass flow rate through the orifice meter was calculated 
from an -equation derived from a hydrogen corrected air calibration of 
the orifice meter. The sample flow rate to the- chromatograph and the 
bypass, flow rate were measured by thermo conductivity mass flow rate 
meters with ranges of 0 to 100 seem (1.39 x 10"^ gm/sec) and 0 to 
3,000 seem (^. 17 x 10"3 gm/sec)-, respectively/ ;The discharge coefficient 
of the injector nozzle, based on orifice meter measurements, normally 
ranged from- 0-73 to 0.7&, with an average value of O.76.. In some 
instances, values of K as low as 0.69 were obtained; calculations 
indicate that a film of dirt on the order of 0.025 millimeter thick 
could have caused this low value. 

Pitot pressure was measured with a 3 .Mi- 7 X 10^ H/m 2 absolute 
pressure transducer and jet total pressure with a 3-^-7 x 10^ N/m 2 
pressure transducer. Survey static pressures were measured -on a 
3.W7 X 10^ N/m 2 absolute pressure transducer. All pressures except 
tunnel-wall static pressures were recorded on automatic balance 
potentiometers. The tunnel -wall static pressures, used to compute 
free-stream Mach number with the known tunnel total pressure, were 
read on mercury monometers and recorded periodically during each 


test run. 



12 


Survey Procedure 

.Data consisting of pitot pressure, static pressure, and volumetric 
hydrogen concentration were taken at downstream locations of 7; 30, 60, 
120, and 200 injector nozzle exit diameters. At each station, data were 
taken at the test, conditions indicated in' the following table. 


Test ; 

condition 

II ■ 
— 3 

30 

6o 

■ 

200 

1 

X 

X 

'X 

, X 


2- 

X ' ' 


’ X 

X 

■* 

3 

X\ ■ 

X 

,x 

X ; 

x 


X 

X 

X 

\ 

X' ’ 

X, 


kb each of the indicated conditions, one vertical and three horizontal 

, t 

surveys were made of the flow field. The vertical survey- was made 
along the jet centerline stepwise from the plate surface outward until 
a zero hydrogen concentration was obtained. Horizontal surveys were 
made at points above the plate corresponding to maximum and half- 
maximum concentration and at a point midway between the plate surface 
and the. point of maximum concentration. 

For each horizontal survey, the positive edge (see coordinate 
system in fig. 1) of the mixing' region was located. A stepwise survey 
was made across the flow field from this point until a hydrogen volume 
fraction of zero was obtained. At each point in the survey a gas 
sample and a pitot pressure measurement were taken. While the gas 
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sample was being analyzed, the probe was moved i to the, next point, and 

the sample line flushed prior to the introduction of a new sample. In 

* 

regions of low pitot pressure, a diaphragm-type vacuum pump was used • 
to withdraw the sample from the tunnel- The sample flow rate was 
throttled to maintain it constant. 

Prior to and periodically during each set of surveys, hydrogen 
was drawn from the supply line and analyzed and the. chromatograph 
readout attenuated to full-scale deflection. Repeatability was checked 
to a variation of less than l/2 of .1 percent full scale. 

Data Reduction and Accuracy 

The raw pressure and concentration data at each survey point of 
each set of surveys - one vertical and three horizontal - were reduced 
to yield values of mass- fraction, 'mixture molecular weight, Mach number,, 
mixture total and static temperature,, velocity, mixture and air mass 
flux, hydrogen -mass flow parameter, and the nondimensional coordinates. 
The molecular weight of the hydrogen-air mixture was computed assuming 
molecular weights of the components of 2.016 and 29*0, respectively. 
Mixture total temperature was obtained from the mixture mass averaged 
total enthalpy computed from the measured total temperatures. The 
mixture was assumed to he- a perfect gas and the values of Mach number, 
static temperature, and velocity computed using the equations for . one- 
dimensional isentropic flow presented in reference 15- Local density 
of the mixture, used to compute the mass flux parameters, was calculated 
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using the ideal gas equation of ’state with the , universal gas constant 
equal to I.986 cal/mole-°K. 

Before each set of surveys was' made, the probe tip was positioned 

^ * « 

along the tunnel centerline on the plate surfaces Vertical surveys ' 

V - , 

were made at this probe location and it was the zero reference point 
for the horizontal surveys which were made by yawing the probe. During 
the course of the tests, it was determined that the centerline of the 
hydrogen-air flow field was not always coincident with the tunnel 
■centerline. This is believed to he a result of the small scale of the 
injector,, prohe tip, and slight asymmetries of the tunnel flow. The 
maximum concentration was therefore sometimes obtained at a point to 
either side of the vertical survey location. When this occurred, the 
y/d location at which all of the horizontal concentration surveys 
"peaked was taken as the true centerline (y/d = 0)- and the vertical 
survey considered to have been made at a point slightly off-center. 

In most eases the distance between the tunnel and flow field center- 
lines, (y/d) a , was less than one injector diameter. Prohe position 
accuracy of the actuator mechanism in the vertical surveys was 
+0. 127 mm, the same as the prohe tip height. This 'corresponds to an 
error in the vertical prohe position of +0.125 injector diameters. 

Since the horizontal surveys were made by yawing the prohe, the 
x/d position of the prohe tip is slightly greater at the edges than 
at the center of the mixing region. At the widest survey location, 
the change in the x position was less than two injector diameters; 



VIII. RESULTS AM). DISCUSSION 
Flow Field Structure _ 

The general structure of the flow field resulting from the normal 
sonic injection of hydrogen in a supersonic air stream is presented in 
figure 6. Flow conditions correspond to a q r of 1.0 which gives a 
value of the jet exit static pressure of 1.1 atmospheres. The how shock 
do™ stream of an x/d of 20 was determined from sehlieren photographs 
and is essentially a Mach line. Its shape is. not 'appreciably affected 
hy a change in q r . Details of the jet structure and "bow shock in the 

i * , ' . , 

vicinity of the injector are not clear because of the small-scale and 
thick boundary layer. Details of the boundary -layer- separation and jet 

J 

shock structure are presented in references- ‘5 , -6, and >8 for boundary 

layers on the order of one injector-diameter thick and, larger Values' of 

* » . » 

» t * r 

q. r than those of the present- investigation. i The severity of- the 

separation depends on the boundary -layer thickness relative., tq the 

1 , * \ ’ 

injector diameter and the- amount- of underexpansion of the jet. • Injec- 
tion through a thick boundary layer, though having- a greater penetration, 

’ i / , i * 

' , * , ‘ 5 

would he turned downstream somewhat before encountering the high-velocity 
mainstream and would result in a weaker bow shock in the free stream. 

The injector hack pressure .would therefore be less for injection- in 
thick boundary layers than, in thin- boundary layers.. Reference 8 used a 
criterion for matched injection by defining an effective back pressure, 
Pfj,- as equal to two-thirds of the free-stream pitot pressure. For the 
conditions in figure 6, p^ is 1.27 atmospheres according to this 
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criterion. Since a sonic injector cannot operate overexpanded, a 
matched pressure condition is the minimum pressure for' operation.- 

Investigation of the injector operation over a range of jet total 

* ' 3 ’ 

pressures indicated that the jet measured. mass flow rate i was -a linear 

- ■ ’ ‘ * ' T 

’ - ' i i ' 

function of the jet total pressure and', therefore,- the jet maintained 

„ i 

sonic operation, down to a value of p, . corresponding to a q - ' of 

* 

approximately Q.b5* This corresponds to a jet exit static pressure, 
and thus an effective back pressure of 0 . 63 atmosphere. For -operation 
in the thick boundary layer of this investigation, with a Mach k free 
stream, the effective back pressure is approximately ko percent of the 
free-stream pitot pressure. It is believed that this result differs 
substantially from that of reference 8 because of the relatively 
thicker boundary layer of the present investigation. 

The data shown in figure 6 are profiles of hydrogen mass fraction 
taken in vertical surveys behind the jet centerline. Trajectories are 
shown for the line of maximum concentration, half -maximum concentration, 
and the point at which the volume fraction is half of 1 percent. As 
can be seen, the hydrogen- jet is rapidly turned downstream try the free 
stream and mixes rapidly in the near field, the maximum concentration 
decreasing to about 12 percent mass fraction in seven injector diameters. 
From an x/d of 7 to x/d of 60, the mixing is slower with the maximum 
concentration decreasing by about- 60 percent. - The major part .of the jet 
penetration into the airstream occurs within seven injector diameters. 
For the thick boundary layer, S/d = 2 . 70 , and low values of q^, the 
mixing region remains almost entirely embedded in the boundary layer. 
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Effect of dynamic- pressure ratio on jet initial penetration .- The 
penetration of a gaseous jet into a supersonic free stream has "been 
greatly discussed in the literature. Zukowski and Spaid (ref. 5) and 
Spaid, et al. (ref. 6 ) define penetration as the point at which maximum 
concentration occurs, while Vranos and Nolan (ref. 7) and Orth, et al. 
(ref. 8 ') consider penetration as the height above the plate at which 
the volume fraction is O.OO5. As used herein, the term, "penetration" 
will refer to the vertical edge of the mixing region and be denoted 
as- (z/d) 0 ; the height at which the concentration- is, maximum is referred 
to as "penetration-to - and will be - denoted (z/d) a . 

A r 

I 

Equations describing the penetration and. trajectories have' 

• * \ t . 

been developed by correlating data .(ref. 7 ) ' are given- here for 
normal, sonic, injection of hydrogen in a, Mach U.03 free stream. 


W 0 - 3-75 q r 0 - 5 (x/d ) 0 - 0866 , 
Wa) a = 3-45 4 r 0 - 535 ( x / d) ; o - 2 5 ? 


( 1 ) 

(' 2 ) 

1 


These equations were derived for data in the far field and are appli- 
cable for x/d greater than 7* Equations (l) and (2) are presented 
in figures 7 and 8 as solid lines with data from the present investi- 
gation. In figure 7 the effect of q r on the initial penetration, at 
an x/d of 7, is illustrated. Equation (2) shows a reasonable agree- 
ment with the data, although it does predict slightly less penetration- 
to -Oftgx at the higher values of q^. A straight line through the 
data points would have a slope of 0.6 compared to 0.533 f° r equation (2). 
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The penetration predicted by equation (l) is about 20' percent low at 
an x/d of 7 and "has a slope of 0.5* A straight line through, the 
penetration data points would have a slope of 0 . 3 . 

Figure 8 shows the. and penetration trajectories for a value 

of q r . of 1.0. Equation (l) predicts an outer boundary of the mixing 
region that increases from about 20 percent low at x/d of 7 to about 
45 percent low at. an x/d of 200. The trajectory predicted by 

equation ( 2 ) approaches the plate at large values of x/d. Data show 
an initial decrease in the a ^g y trajectory to an x/d of 30 but an 
increase at stations farther downstream. Equations (l) and (2) were 
obtained for helium injected from a 0 . 478 -centimeter -exit diameter 
nozzle mounted normal to the wall of an 8 . 89 -centimeter-diameter duct. 
The curvature of the wall and the possibility of interference from the 
opposite side of the' duct could account for the discrepancies between 
the present data and the equations. Another correlation for the pene- 
tration of a normal jet and reported in reference 9 is 

(z/d) = 1.446 q 0*392(M/M.) 0 - 6l 3(x/d) 0 -°39 6f0 ;3T9;M j /3^ ( 3) 

For sonic injection in a Mach 4.03 free stream, the equation reduces 
to 

(^)„ = 3A0 13,1 ’ ft) 

Equation ( 3 ) was derived from data taken at free-stream Mach numbers" 

J- . 

of 1.6 and 3-0 over an x/d range of l4 ,to 167 . The- injected gas was 

. ■ 1 ( i » 

* 4 \ 

simulated methane. The effect of boundary -layer thickness was not • 
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investigated. Equation (4 ) is presented in figures 7 and 8 as dashed 

\ > 

lines and ’gives a tetter prediction of the' penetration than equation (l). 

\ 

An equation derived from a correlation of . the present data at values of 


x/d less than 120 - shown as a -dashed line in figure 9’( c ) ** is 


/ 0.300 ’ :0‘.l43 

(z/d.) 0 = 3*87 <1£\ ; (x/d) 


(5) 


The slight- differences in the constants of equations (t),and (3 ) • could 
he due to the different injected gases and different' boundary layers. 

t ‘ * • * 

The broken lines in figures 7 and 8 are the results from 'a method 
presented in reference 8 for calculating the jet centerline trajectory 
(“max trajectory) in the near field. This method considers the 
emerging jet as being composed of cylindrical elements of length d(h/d) 
with the aerodynamic drag on each element computed from empirical 
equations. The equation from reference 8 is 


6.91 



ae 

c D (e) sin2(e) 




( 6 ) 


where 


Cp(0) = 1.2 + (Mqo sin B)^^, 0- < Moo sin 0 < 1 (7) 

C D (e) = 1.06 + l.l^Moo sin Mco sin 8 >1 (8) 

The integral is evaluated hy selecting a value of © and integrating 
to the initial value, 0^. For normal injection, 0^ is 9°°* 
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Equation (6) underpredicts the effect of q r on the initial pene- 

i 

tration at an x/d of 7, "but again, the effect of boundary layer is 

. * v- 

not accounted for. Calculations were made using the mass averaged 
boundary-layer conditions rather than free stream in equation* (’6) with 

t * - » ■> 

no significant change in the trajectory r or the effect of q . * 

Effect of dynamic pressure ratio on the penetration trajectories . - 

Figure 9 presents the trajectories of maximum concentration, 1 half- 

» * 

maximum concentration, and the penetration height, correlated with the 

jet-to-free-stream dynamic pressure ratio; For all -/values' of q r , 

except a value of 0.5, the penetration-to - a. _ decreased with increas- 

' * ’ < 

ing x/d to an x/d of about beyond which penetration-to - 
increased. The minimum value of (z/d) a occurred farther downstream 
the higher the value of q r - In order to bring these minimum points 
together a factor q r ~ ‘ was applied to the x/d coordinate producing 
the family of curves presented in figure 9(a). The fact that the maxi- 
mum concentration trajectory turned beyond parallel with the plate 
surface, as evidenced by the initial decrease in (z/d ) a , is thought to 
be caused by the thick boundary layer. Figure 9(b) presents the tra- 
jectories of the point at -which the mass concentration is one-half the 
maximum at each x/d station. At all downstream stations the (z/d)^ 
coordinate was taken in that part of the flow field between the point 
of maximum concentration and the edge of the mixing region. The effect 
of q^ on the coordinates of these points is such that a factor. of 
c r " u ‘ applied to the x/d coordinate produces a family of curves 
similar to those for the (z/d) trajectories. The family of curves 
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for the (z/d) a and (z/d)^ trajectories "both show two distinct 
characteristics. For Ojj^ trajectory family , the data points essen- 
tially lie on a straight line of negative slope at values of the x 
coordinate less than 30. Downstream of this minimum point the 
trajectories diverge with the slope increasing for increasing values 
of q_ T - The half -maximum trajectories exhibit similar properties with 
the dividing point, for a q r of 1.0, occurring near a value of x/d 
of 15* At all x/d stations, an increase in q r produced a propor- 
tional increase in (z/d) Q such that the individual trajectories were 
brought together by multiplying the (z/d) D coordinate by q r ^ as 
shorn in figure 9(c). The solid line is a fairing through the data 
points^ the dashed line is a straight-line approximation to the pene- 
tration trajectory at values of x/d less than 120, and is represented 
by equation' ( 5 ). 


Decay of maximum concentration . - The decay of the maximum concen- 
tration as a function of x/d-j and q r is presented in figure 10. 

The effect of -,q was, found to be such that a factor of q 
1 r 

multiplied by the x/d-. coordinate provided a reasonable correlation. 
Downstream of ^an x/d of 30 the concentration decay may be repre- 

t , ' * * 

• 'Rented by ,a straight line with a slope of -0-8; that is, the decay of 

V f ’ Q g 

is inversely, proportional to (x/d^-) ’ . Extrapolating this 

j f 

straight line - ' shown' as a dashed line in figure 10 - to a value of 
" a max ’ yields a length equivalent- to the potential core in 

. ^ i ’ * , } 

coaxial flow. For values of q_^ from 0.5 to 1 . 5 , the equivalent 
potential core length ranges from 1.2 to 2.2 jet diameters. Compared 
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to coaxial flow the mixing for normal injection- is much faster in the 
near field, yet slower in the far field. As reported in reference 2, 
the concentration decay in the far field for' coaxial mixing i-s inversely 
proportional to x 2 . 

Presentation of Data 

Profiles of concentration, velocity, and total pressure were non- 
dimensionalized for the vertical survey and the horizontal survey 
through the point of maximum concentration at values of q ■ of 0-5, 

1-0, and 1.5 and all x/d stations. These profiles are presented in 
figures 11 through 1J> for the vertical surveys and figures lb through 
16 for the horizontal surveys. 

Vertical survey profiles . - For each vertical survey the, largest 
value of mass concentration, arid the corresponding value of z/d 

were obtained by fairing a curve through the data -points. The value 
of d is not necessarily the maximum concentration due to the uncer- 
tainty of aligning the sample probe with the jet centerline (see the 

Data -Reduction and Accuracy Section). 

* • ' 

The- [Concentration profiles normalized by a re f and with the 

t 

origin of the coordinate system shifted to (z/d) = (z/d) re £ are P re “ 
seated -in figure 11 for each value of x/d. At values of z/d greater 
than the value of z/d at <x re £, the vertical coordinate was non- 
dimensionalized by the distance from the point of a ^ to the edge 
of the mixing region, (z/d) Q - (z/d) At values of z/d less than 

, (z/d ) £, the vertical coordinate was nondimens ionaliz ed by (z/d) ref . 



23 


For z/d greater than- (z/d) rg £, a change in the value of q r has 

negligible effect on the shape of the profile. However, the profile 

' * 

.shape does change with x/d. To .compare the shape of the profiles at 

* ‘v 

values of z/d. greater than ( z/d) re ^ as the mixing progresses down- 

* , . \ 

stream,- a Gaussian-type exponential curve of the form 


a 


a ref 


=- exp / -5 


z/d - (z/d) ref 


(z/d) Q - (z/d) 


ref 


( 9 ) 


is .shown. for each x/d. The value of b-j,, the vertical profile shape 
index, required to give a reasonable fit of equation (9) with the data 
is given in the tables in figure 11 . At x/d stations downstream, of 
30, the value of is constant at a- value of. 2.70, indicating that 
the flow field has become fully developed. This may be associated with 
the fact that the maximum, concentration trajectory in figure 9 has a 
minimum near a value of (x/d)q r ' of 30 and that the maximum con- 
centration decay in figure 10 has a constant slope downstream of an 
x/d of 30. 

The portion of the profiles at values of z/d less than (z/d) re f 
show no predictable effect of q r . However, at values -of x/d less than 
60 the slope of the lower half of the profiles increases, in general, 
with an increase in q^. That is, in the vicinity of the plate surface 
the mixing region is generally more uniform at the lower values of q r > 
The slope of the profiles near the plate surface decreases with increas- 
ing x/d. At an x/d of 60 and above, curves faired through the data 
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points would intersect with the plate surface at values of a/a re £> 
between 0.6 and 0.7 for all values of q_ r . 

The nondimen 3 ional velocity profiles are presented in figure 12 
for all x/d stations at each value of q r - Also shown for comparison 
is the undisturbed "boundary -layer profile taken at the injector exit 
station (see fig. 2). For each data profile , the value V Q is the 
velocity at the edge, of the mixing region at a height (z/d) Q above 
the plate. For the boundary-layer profiles V Q = Y m and (z/d) 0 = 6. 
At x/d stations of 7 and 30, the effect of increasing q r can be 
seen as an increase in the peak velocity located within the profiles. 
This is due to the fact that for constant initial conditions the mass 
flow rate of injected gas increases directly with q. r so that there - 
is more high-velocity hydrogen near the injector at the higher values 
of q r . The initial jet-to-free-stream velocity ratio is approximately 
2. At values of x/d downstream- of 30, the shape of the velocity pro- 
files approach that of the boundary layer, the higher values of q r 
-requiring a greater length. 

*4 * % 

'The quantity , z-^ . in the tables of figures 12 and 13 is the height 
of the undisturbed streamtube upstream of the inj ector that contains - 

l 

the same mass flow of air as the mixing region downstream. It was com- 
puted by assuming that the undisturbed streamtube has the same cross- 

i ' , 

, 'sectional shape as the mixing- region and will be discussed later in 

a v 1 , * 

'more detail! For the boundary -layer survey, z-^ = 6. Each profile 
-then represents a segment of boundary -layer-type flow with a thickness 
of . z-^ when referenced to the undisturbed flow. 
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Nondimen sional total pressure profiles corresponding to the same 
conditions as the velocity profiles are presented in figure 1J. The 
boundary-layer total pressure profile from figure 2 is also included 
for comparison., As mth the velocity profiles, the total pressure 
profiles at stations far downstream approach the same shape as the 
boundary -layer profiles'. However, near the edge of the mixing- region, 

i 

at large' values of x/d, the slope ef the total pressure profiles is 
larger than the boundary -layer profile. The total pressure at the 

edge of the mixing region generally increases with increasing x/d. 

« * * f 

-This \is due to. an addition to the mixing region of free-stream air 
that has .passed through a weaker portion of the how wave. All of the 
profiles exhibit the same degree of total pressure "loss due to the 
loss in momentum required to turn and accelerate the injected hydrogen. 
This region of low total pressure extends over approximately 60 percent 
of the vertical height of the flow fields and is most severe at the 
upstream stations. 

Horizontal survey profiles .- Nondimensional concentration profiles 
for the horizontal survey through the point of maximum concentration - 
are presented in figure it. Here the concentration is normalized by 
the maximum value which occurred a distance (y /d ) a from the tunnel 
centerline, or the point at which the vertical survey was taken. The 
lateral coordinate* is nondimensionalized by the average distance from 
the centerline to each of the lateral edges of the mixing region, 

(y/d) Q . The values of (y/d) 0 have a somewhat random variation . 
varying from roughly 5 to 8 from an x/d of 7 to 200 with a t20 percent 
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deviation at, a .given station. As with the vertical surveys, the dynamic 

4 * 

pressure ratio has little or no effect on the profile shape at each x/d. 
To compare the effect of x/d ■ on the profile shape, a Gaussian-type 
exponential curve is presented. The equation is of' the form 


a ■ 
c^max 


= exp 



y/& - (y/ d ) a 

r^r. 



( 10 ) 


where bg is the horizontal profile .shape index. The values of bg 
required to give a fit of equation (10) with the data are given in the 
tables in figure lit- . At values of x/d of J , 30, and 6o, bg is 
constant at a value of 1*5 and increases to a value of 2.0 at an x/d 
of 200. Equation .(l0) fits the data reasonably well, especially at an 
x/d of 120 and 200. 


The velocity profiles for the horizontal surveys are presented in 
figure 15. The velocity, V Q , used to nondimensionalize the data points 
has the same value as that used for the vertical surveys in figure 12. 


For absolute values of 


[y/ d - (yAOq] 
(y / d ) 0 


greater than 0.5, the profiles 


are fairly uniform with the level' of the velocity ratio, more a function 
of the value of z/d at which the survey was made than it is x/d. 
Generally, in the outer part of the flow, field the velocity increases 
with increasing z/d. I-Tear the center of the mixing region - absolute 
'values of the y coordinate less than 0,5 - the velocity level 
decreases with increasing x/d. Thus, as the mixing progresses down- 
stream, the velocity profiles, instead of having a maximum near the 
center of the flow field, have a minimum. This is a result of the 
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flow field "becoming established at large values "of x/d as indicated 
in figure 12 "by the vertical velocity surveys approaching the- shape 
of the undisturbed boundary layer. 

Profiles of total pressure for the horizontal surveys are presented 
in figure 16. As with the velocity profiles, the total pressure pro- 
files exhibit different trends near the edge of the flow field than 
near the center. Wear the edge of the mixing region, the level of 
pressure increases with increasing z/d, while in the vicinity of the 
centerline, the pressure level is a function of x/d and, in general, 
increases monotonically with x/d. As the mixing -progresses downstream,' 
the -pressure profiles become steeper as a result of higher momentum air 
being added to the mixing regions at the edges. The; pressure level at 

the centerline never recovers from the injection disturbance and is 
. « 

‘less than 8 percent of free-stream total pressure even at 200 injector 
diameters' downstream. 


Flow Field Contours , 

$ Contours of mass fraction, hydrogen mass flow rate parameter, and 
ai-r mass flow parameter were obtained by cross -plotting the vertical 

r 

and horizontal prpf lies at constant values of -the particular parameter. 

* > - ' ' 

Contours^ of a and 3 are presented in figures 18 and 19, respec- 

’ * * m 

tivelyj contours of' 6 are not presented, but the results of the inte- 
gration of the contours performed as a continuity check and an indica- 
’■tion of the overall accuracy of the data are given in figure 17 • The 
contours of fuel mass flow were plotted in the nondimensional form. 
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fi/6 max where is the maximum value for each set of surveys and 

t 

generally occurred 'near the point of g ^y . In terms of £, the inte- 

S , ' 

grated -mass . flow rate of hydrogen is 


m. , 

xnt 


^ 'C p v)^ max 


r 

J o 




( 11 ) 


Dividing equation '(ll) ty A c 


raitej 


m. = 
3 


rt (pVWd,2 


gives 


and the measured hydrogen mass flow 


= >Wo r l '° ( i \J a\ 

mj itdj 2 J o \^max/ vW 


( 12 ) 


as the ratio of integrated-to-measured hydrogen mass flow. Results of 
evaluating equation (12.) at the various downstream stations are .given 
in figure 17. The solid symbols in figure 17 represent the average 
deviation of equation (l 2 ) from unity. As can be seen., the accuracy . 
of the data is better at the far downstream stations. At all values 
of x/d, the deviation of equation (12) from unity increases with 
increasing values of q r . This is probably due to the larger gradients 
in the concentration and velocity - which are associated with the local 
turbulence level - produced by the stronger jet at the higher values 
of q r . The characteristics of a binary gas flow field that affect 
the accuracy of concentration measurement are discussed in reference 3. 
Accuracy of concentration measurements in a flow field composed of 
■gases with a large difference in molecular weights- is also affected- 
by the sampling prohe-tip geometry. An investigation reported in 
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» i. 

reference 14 for subsonic flow of coaxial jets indicated that a probe 
tip with a rapid internal expansion provided better results than one 

T * * 

with a -long length of constant internal area when compared to measured 

L 4 I , ' 

mass flow rates. Reference 3 indicated that the probe and sampling 

lines should be' free of any obstruction so that the probe tip will 

f , i , 

operate at nearly full capture. For flow fields of this nature, 
■differences between- integrated and measured mass flow rates of 20 per- 
cent are considered typical. 

The contours represent cross sections of the flow field in the 
y-z plane and are hounded by the contour for a concentration of zero 
which defines the outer edge of the mixing region. Above a value of 
z/d of about 3 the zero concentration contour may be approximated by 
a semicircle with the center at the point of maximum concentration. 
Below a z/d of approximately 3 the zero concentration contour rapidly 
spreads laterally in the boundary layer. 

For each value of q r the overall width of the a = 0 contours 
is essentially constant at x/d stations downstream of 30. The height 
of the contour, which is approximately the penetration, almost doubles 
from an x/d of 7 .to an x/d of 200. The contours for a particular 
value of a other than a = 0 decrease in area with increasing x/d 
and decreasing value of q r as shown in figure 18. 

The air mass flow rate contours in the nondimensional form P/P max 
and- contained within the a, = 0 contour are presented in figure 19* 

The intersections of the P/P^ax coir k our and the a = 0 contour were 
determined by plotting the value of 0 at the point in the horizontal 
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surveys where a '= 0 as a function of z/d and interpolating to find 
the coordinate corresponding to a particular value of p/ fWx* ® ie 
mass flow rate of air contained within the mixing region was obtained 
by evaluating the integral 


ol.O ‘ 

“1 = IW*o (s/fJm axWW 

J 0 


(13) 


where .A is the cross-section area bounded by each P/3 ma}C contour 
and A q is' the cross-section area within the a = 0 contour. Results 
of' the integrations were used to determine the average fuel-air ratio 
and the size of the undisturbed streamtube upstream of the injector 

t f 

that contains the same air mass flow rate as the mixing region 
downstream. ’ 

" The undisturbed -streamtube area,, A-^, was . calculated for -each 
x/d and q by assuming that it was of the same shape as A Q ’. That 

is., the .ratio of the height, z-^,, to the average width, y^ ■= A^./z-^, 

* v , ' 

of the streamtube was assumed constant. This gives the following 

' > i 

t * 

i equation for A^: 

x2 


\ = wv 


(a) 


The height of the undisturbed streamtube was calculated from the con- 
tinuity equation and the definition of boundary -layer displacement - 
thickness 

m-L = [ (pV)dA = (pV^Czq " S*)^ 


(15 )• 
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but 


Tl = kjz ± = A o z x !z o 2 (16) 

’ i 

Substituting equatipn (l6).into (15) and solving the resulting quadratic 
for z x gives * 




(17) 


Values of Zq obtained from equation (17) for the integrated air mass 
flows of equation ,(15) are given in the tables of figures 12 and 15. 

In .application to the design of supersonic combustors, the size and 
shape of the undisturbed streamtube provides information as to the 
spaaing of injectors and the injector size to obtain penetration to the 
combustor centerline. If a stoichiometric value of f is desired the 
combustor entrance must have a half -height equal to the value of z x 
that corresponds to the x/d station at which a stoichiometric average 
fuel-air ratio was obtained. Confining the flow field by the addition 
of an opposite wall would he expected to change the mixing rate and 
penetration from that presented herein hut. would yield the same average 
fuel -air ratio. The injector spacing is the average width of the 
undisturbed streamtube, y^* Ah estimation of these parameters is dis- 
cussed in the next section. 
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Estimation' of Cold glow Mixing Parameters 
'The mixing length relative to the combustor entrance height and 
the injector spacing required to give a stoichiometric average fuel-air 
ratio (f = 0.0293 for hydrogen in air) has been estimated by super- 
imposing the flow field cross section of a single injector. A schematic 
for a two-dimensional configuration with staggered injection from both 
walls is presented in figure 20. It is assumed that a particular con- 
centration contour from the single jet flow field coincides with the 
same contour of the adjacent and opposite injectors. For the purposes 
of this discussion, it is assumed that the superimposed flow field is 
coincident along the contours of half -maximum concentration, shown as 
dashed lines in figure 20. It follows that the air mass flow that 
mixes with the hydrogen from one injector is contained within the 
half contour, since the area hounded by the adjacent contours 

covers the entire cross section of the combustor. Denoting the area 
contained by the half contour as A^, equation (ij), for the air 

mass flow, becomes 


■h = fma ^ 5 / 0 <e/<W> d (VV ( 18 ) 

The height and width of the confined air streamtube are given by 
equations (17) and (l6 ) by replacing A Q and z Q with Ae> and z^, 
respectively. ‘ 

Values of average fuel-air ratio 'obtained from' the total injected 
hydrogen mass flow and the integrated air mass flow for the unconfined 
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jet (eq. (13)) and. the 0'/cc niax of 0.5 contour (eq. (l8)) are presented 

in figure 21 as a function of x/d-j and q r . The effect of q_ r on 

the decay of f for both conditions was found to be such that a factor 
-l/2 

of q y 1 provided a reasonable correlation. For the unconfined jet 
the value of f Q is below stoichiometric for all x/d stations con- 
sidered. Downstream of 30 injector diameters, the decay of f Q is 

. , . 0 . 6 - 

inversely proportional to approximately (x/dj ) . . For the case of 
simulated opposite wall injection the average fuel-air- ratio, f , 

i 2 

is considerably higher than f Q and becomes stoichiometric at an 

x/d. of 200 for a a of 1.0. It is apparent that the lovest value 
1 r t . . * 

of q r will provide a- stoichiometric average fuel-air ratio in the 

shortest distance. Also, using a coincident concentration contour with 

' * 7 ■ ~ 

a value less than half maximum would, yield 1 a stoichiometric value of ■ 

~ , « , * 

f at a shorter x/d since the resulting curve for the decay of f 

would lie somewhere between the curves for f Q ' and ’ f c* However, the 

* , 

uniformity of the flow field would need to. he investigated to select 

the optimum value. 

To determine the combustor size, the average fuel-air ratio of 
the simulated combustor flow field is plotted in figure 22 as a func- 
tion of x/z-^. The value of represents the combustor entrance 

half -height. From figure 22, f^ is stoichiometric at a value of 
x/zj- of approximately 40. The hydraulic diameter of a two-dimensional 
duct is twice the duct height, or 4z^. The mixing length required for 
a stoichiometric average fuel-air- ratio is approximately 10 hydraulic 
diameters, which agrees with values used for engine design. For a q r 
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of 0-5, the value of for a stoichiometric value of f~ is 

approximately 3.3 injector diameters - corresponding to an x/d- of 
1^-0 from figure 21 . The value of y^, -which is the required injector 
spacing corresponding to these conditions, is approximately 3*0 injector 
diameters. 



IX. CONCLUDING REMARKS 


Comparisons of the present data to data correlations from, other 

sources indicated a significant effect of the thick boundary layer on 

■the initial, penetration of the jet, the downstream trajectories, the 

jet effective back pressure, and probably the mixing rate. All the 

data correlations underpredict'ed the effect -of jet -to -free -stream 

dynamic pressure ratio on the initial penetration and the penetration 

trajectory. Equation (3) from reference 9 provided the best agreement 

0.3 

to the present data which correlated as a function pf q_ r . The 
trajectory for the penetration -.to - was initially turned beyond 

parallel to the plate- surface and then increased monotonic ally with 
the slope of the trajectories increasing with dynamic pressure ratio. 
Minimum values of (z’/d) a ‘.occurred at values of x/d ranging from 9 
to 45 for values of q r between 0..5 and'1.5, respectively.. As a 
consequence of the overturning of the maximum and half -maximum concen- 
tration trajectories, simple empirical data correlations are not 
possible. 

Investigation of the jet operation over a range of pressures 
indicated that the jet effective back pressure - the minimum jet exit 
static pressure for a sonic injector - for operation in the thick 
boundary layer of this investigation is approximately t-0 percent of 
free-stream pitot pressure. 

From the level of the measured maximum concentration near the 
injection station, it can be concluded that the turning and initial 
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mixing of the hydrogen jet is very rapid. Downstream of an. x/d of 30, 
the decay of the maximum concentration is' inversely proportional to 

o 3 » 

(x/d_j ) * , which is a slower mixing rate than coaxial mixing in the 
far field. Mixing, in the near field is much faster for normal injection. 
The equivalent, potential core lengths for normal injection were esti- 
mated to range from 1.2 to 2.2 jet diameters for values of q r between 
0-5 and 1.5, respectively. 

Examination of concentration profiles oh the vertical centerline 

, * T 

1 ’ 

suggest that the profile shape, above the point ”of maximum, concentration 
is not dependent on the dynamic pressure ‘ ratio, and may be represented 
by a Gaussian -type function. The profiles show similarity at values 
of x/d equal to or greater 'than 60 and are nonsimilar- at values of 

t * 

x/d of 7 and 30* Horizontal concentration profiles through the point 

' , . * 
of maximum concentration are also represented' by a ,Gaus sian-type 

' ' , . 1 1 

function -and exhibit similarity at value's of' x/d- -less than 60. 

■1* 

The loss in momentum and -total pressure of the airstream, caused 

! 1 

t 

by turning and accelerating the hydrogen jet, resulted in an extensive 
region of very low total pressure. The region extended over 60 percent 
of the height near the center of the flow field and 4-0 percent of the 
width with a total pressure less than 10- percent of the free stream 
even at 200 injector diameters downstream. 

Investigations of the mixing patterns at each, station to obtain 
overall data discrepancies as judged' by the measured fuel flow indicated 
a general trend of large discrepancies at upstream stations where values 



3 ? 


of the maximum concentration were high, and discrepancies approaching 
zero at downstream, stations where was near stoichiometric. 
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Figure 1.- Sketch, of model. All dimensions, in centimeters 



Figure 2.- Boundary layer surveys on flat plate model at injector exit 

station $ no injection- 
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Figure 3.- Tunnel and secondary flow schematic 
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□ 


a max 
v « 0.'005 
Eqs. 1 & 2 
Eq. 4 
Eq. 6 


Figure 7*- Effect of dynamic pressure ratio on the jet penetration , 
x/d = 7j comparison of., data * and correlations. 



Figure 8.- Maximum concentration and penetration trajectories 
= 1.0; comparison of data and correlations. 
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(a) Maximum concentration trajectory. 



(x/d)q r ' 8 


(t>) Half -maximum concentration trajectory. 



(c) Penetration trajectory, v = 0.005* 

Figure 9 .- Effect of dynamic pressure ratio on penetration trajectories. 





(a) x/d = 7, 30, 60. 

. Figure 11.- Mbiidimensional concentration profiles, vertical survey. 
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Figure 12.- Nondimensional velocity profiles, vertical survey 
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Figure 12.- Concluded 
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Figure 13- ~ 


No ndimen sional total- pressure profiles, 


vertical survey. 
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Figure 13- - Concluded 



(a) x/d =3 'J } . 30., 60. 

Figure lU.- Hondimen sional concentration profile s, horizontal survey 
through point of maximum concentration. 
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Figure lk.- Concluded. 
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l7T% 


(a). q. r = Oi'5* 
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(b) q. r = l.o. 

Figure 1 6 .- Efondimensional "total pressure .profiles^ horizontal survey 
through, point of maximum concentration. 
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Figure l6.- Concluded. 
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(t) q. r ® 1.0 (concluded). 


Figure 19 •- Continued. 
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Jc) a r = 1.5. • 

Figure 19 Continued., 












Figure 20.- Schematic of flow field cross section for opposite staggered 
injection by superposition of single jet contours- 
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Figure 21.- Average fuel-air ratio decay. 



Figure 22.- Decay of average fuel-air ratio for simulated opposite 

wall injection. 



